CLNS 97/1472; 



hcp-th/9703212; 



A ID Model for N-Level Atoms Coupled to an EM Field 



Zorawar S. Bassi 
Andre LeClair 

Newman Laboratory 
Cornell University 
Ithaca, NY 14853 



We construct a model for n-level atoms coupled to quantized electromagnetic fields in a fibrillar geometry. 
In the slowly varying envelope and rotating wave approximations, the equations of motion are shown to satisfy 
a zero curvature representation, implying integrability of the quantum system. 
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Introduction 



The interaction of radiation with two-level atoms has been extensively studied under various approxi- 
mations. In one spacial dimension, the reduced Maxwell-Bloch equations resulting from the slowly varying 
envelope and rotating wave approximations are known to be quantum integrable [1]. In this paper we 
generalize the one-dimensional case of two-level atoms to that of n-level atoms. 

In the first of two parts, we construct the fully quantum n-lcvcl model. The system consists of n-level 
atoms distributed in a fibrillar geometry, interacting with radiation through a minimally coupled hamiltonian. 
In the remaining section we apply the approximations, and show that the Heisenberg equations of motion 
for the reduced system satisfy the so-called zero curvature representation. This implies that the system is 
integrable and can be solved by the inverse scattering method. 

0. Mathematical Background 

Let us first recall the sl n Lie algebra. The n 2 — 1 generators, written as 

{Eij,H m 1 1 < i ^ j < n, 1 < m < r}, (0.1) 
where r = n — 1 is the rank, satisfy the following brackets (in the Chevalley basis) 



[Eij , E, 



kl 



SkjEu — SuEkj, if SaSkj = ' 
j'-i 

^2 H m (i < j), if SuSkj = 1 



(0.2a) 



[E^, H rn j — (5j m — Si m — 5j m+1 + <5j m+1 ) E^ , (0.2b) 
[H a ,Ht]=0. (0.2c) 

The set spanned by {H m } being the Cartan subalgebra. (Note that in the bracket (0.2a) a term of the 
form Si^jEkk is formally equal to zero, even though Ekk has not been defined. The set {E^ \ 1 < i, j < n} 
satisfying the first relation in (0.2a) is a basis for the algebra gl n .) A representation p of sl n will be denoted 
as 

{Ee j =p(E ij ),H^=p(H m )}. (0.3) 
The r x r Cartan matrix A has the explicit form 

A U v — 2<5 Mt , — ^utj-i ~ ^md+i- (0-4) 

It is a symmetric matrix with diagonal elements 2 and nearest off-diagonal elements —1. We now proceed 
to build our quantum system. 



1. The Interacting N-Level Hamiltonian 

We model a free n-level atom as having a single electron with eigenstates |i), i = 1, 2, . . . , n, and energies 
ei > £2 > • • • > £«• The energy splitting between states |i) and |i + 1) will be denoted by Wj or u)u +1 (h = 1) 

Wi = Wjj+i = ej - ej+i, 1 < i < n. (1.1) 

The second notation can be generalized as follows 

uiij=ei-ej, l<i<j<n. (1.2) 

The notation Wj is only defined for the energy splitting between successive states. 

To define various atomic operators, we first introduce fermion creation and destruction operators {bi, b\} 
for 1 < i < n. The operator b\ (bi) creates (destroys) an electron in the i-th level. These operators satisfy 
the til^cbrEi 

{b i ,b j } = {bl,b]} = 0, {b i ,b]} = S ij . (1.3) 
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The atomic operators can now be written as 

Oij=b\b h l<t,j<n (1.4) 
or linear combinations of the O^s. The action of Oij on an atomic state |k) is given by 

O ij \k) = b\b j \k)=6 kj \i). (1.5) 

From (1.3) the general commutator for the O operators is 

[Oij,O k i] = OuSjk - OkjSu. (1.6) 

Operators of the form Oi<j are referred to as raising operators. These cause a transition from the lower 
energy state |j) to the higher energy state |i). Similarly the operators Oj>i are lowering operators. We also 
define a set of commuting operators, denoted H m , as follows 

U m = O mm - O m+lm+1 , 1 < m < r. (1.7) 

The set {Oi^j,Ti. m } satisfies (0.2), thus forming a representation of sl n . (Note: We shall often use the 
notation X a< b (X a> b) to mean X a b with a < b (a > b) for X any quantity, operator, c-number, etc.) 
By appropriately choosing the arbitrary lowest state energy e n to be 

r 

Zn = ~ A rrn^m, (1-8) 

rn—1 

where A~ x is the inverse Cartan matrix, the free atomic hamiltonian can be written as 

H ,tom = J- A- U IUJ V U U . (1.9) 

l<u,v<r 

To couple the atom to an electromagnetic field we make use of the minimal coupling prescription. The 
standard hamiltonian is 

H = H+ + ff atom + Hint, (1.10) 

where 

Hr m = ^- e P-p+v(x) (i.ii) 

Hint = "2^(P ' A W + A*) ■ P) + ■ ^(*)> (1-12) 

and Hq is the free field hamiltonian. If the spatial variation of the vector potential A is small across the 
atom, we can take its value at a fixed point x inside the atom. Using 

P = ^[x,Ht m ], (1.13) 



we get 



A*) + A*) ■ P\ h ) - - e a )A(x ) ■ (a|d|b), (1.14) 



where d — ex is the electric dipole operator. Since d is a vector operator and the atomic states are assumed 
to be parity eigenstates, we have = 0. For i < j, the matrix elements are of the form 

(i|d|j) = d i:i <: in '>n. (j\d\i) = d i:j <: in '>n. (1.15) 
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where dij > and the unit vector n gives the spatial orientation of the atom. This shows that the dipolc 
operator can be expanded in terms of the raising and lowering operators as 

d = nY J {dije iaii O ij + rfyft u ^O :ji ). (1.16) 

i<j 

The interaction hamiltonian becomes 

H int = -iA(x ) ■ n^wydtf (0«e ,a « ~ 0jie _,aw ) + - — A(x Q ) ■ A(x ). (1.17) 

To reduce this system to a one-dimensional model, we make use of the fibrillar geometry. The atom 
can be thought of as an impurity in an optical fiber of cross-sectional area A and length L, with L » \[A. 
Taking the fiber along the x direction, the reduced field action is found to be (see [2] for details) 

1 f 1 

^ ^Maxwell = dxdt -(dt<pd t <i) - d x <j)d x (f>) , (1.18) 
where is a dimesionlcss scalar field defined through 



I-S = V^. (1.19) 

Here A is the vector potential depending only on the x coordinate and A c a is the effective fiber cross-sectional 
area. The field <j> satisfies the commutation relation 

[<t>[x, t), d t <j>{x', t)] = i5{x - x). (1.20) 

From the action the free field hamiltonian is found to be 

Ht = I dx \ [(d t cb) 2 + [d^f] + -^£-0 2 (xo), (1.21) 

where the last term is the quadratic potential term taken from H mt . Now the interaction Hamiltonian is 

ffint = -^(so)5>y/M0«e ia " ~ O jie - ia «), (1.22) 

i<j 

where (explicitly showing % and c) 



i. ,1^1,. (1.23a) 



For 1 < m < r we define 



Pm — ^ c j^ — d mm+ i — Pmm+l- (1.23b) 

The (3 parameters are the important dimensionless coupling constants of the model. The spontaneous decay 
rate T\- — 1/t,?- of a single excited atom from the state |i) to the state |j) is given by 

I« = ^W (1-24) 

Next to make the transition to a continuous system. For N atoms positioned at x = x m , m = 1, . . . , N, 
let d m — e(x — x m ) and Oij(x m ) be the dipole and transition operators on the atom at x m . The matrix 
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elements of d m are independent of the position, however, the orientation can vary from atom to atom. The 
operator d m can be written in terms of the single atom matrix elements as 

d m = Sm^iiij^We^ + O ji (x m )e- iai i). (1.25) 

i<j 

For simplification, we consider the situation where all atoms are aligned n m — n (e.g., by an external electric 
field), giving 

#int = -\f dx<t>{x)Y,UijPiA ii(x,tY aii - 3l {x,t)e' 1 ^), (1.26) 
J i<j 

where we have introduced the space-time dependent transition operators 

N 

t)S{x-x m ). (1.27) 

m— 1 

The discrete operator Oij(x m ,t) acts only on the atom at x m to cause a transition from |j) to |i). The 
Ti m (x,t) operators are defined similarly. The general commutator for the space-time transition operators is 

[O ij (x,t),O k i(x',t)] = (Ou(x,t)S jk -O kj (x,t)S u )S(x-x'), (1.28) 

from which it is easily seen that the algebra (or now more appropriately current algebra) satisfied by the 
set {Oi^j(x,t),H m (x,t)} is identical, aside from the delta function factor, to the sl n algebra (0.2). The free 
atomic hamiltonian takes the form 

H* tom = Jdx K^vH u {x,t). (1.29) 

l<u,v<r 

The complete hamiltonian for the system is therefore H atom + H$ + H lnt , with i?J tom , H$ and H int 
given by (1.29), (1.21) and (1.26) respectively. 

2. Two Approximations and Integrability 

We now make use of two approximations common in quantum optics to further simplify H mt . These 
being the slowly varying envelope and rotating wave approximations. 

In the slowly varying envelope approximation, we assume that near resonant photons with energies 
w uoij are most relevant. Then the scalar field </> can be expanded about the various resonances as 

<f>(x,t) « ^(e-^'-^OM) +e i "^ t - x %(x,tj), (2.1) 

i<j 

where ipij{x, t) and ipjj{x, t) are destruction and creation fields with mode expansions 

1 f dk 

^(x.O) = -= J -^a u (*«)e*-, (2.2a) 

and 

a i3 (k e ) = a lj {k e + ojij), alj(k e ) = alj(k e + w„). (2.3) 

Here (at-) is the usual photon destruction (creation) operator. (Note that we are only considering 
right-moving plane waves.) The operator a,ij(k e ) (a\j(k e )) destroys (creates) a photon with energy 

\k\ = | fee + \ « Uij. (2.4) 
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Thus k e acts as an "envelope" vector about the resonance. The photon operators satisfy the standard 
commutator 

[ay(fc),ai I (fc / )] = [atf(*0>4i(*')] = fcfcfy*(* - *0- (2-5) 

All other commutators vanish. In writing (2.1) and the above commutaion relations, we have assumed that 
all resonances w y {j{n 2 — n) in total) are distinct (i.e., uJij = uiki i = I and j = k) and sharp. From 

(2.5) the component fields satisfy 

[ipij(x,t), ^ kl {x',t)\ = -^—Si k Sji5(x - x'), (2.6) 

ZU>ij 

with all other commutators zero. 

The rotating wave approximation reduces the number of interactions in H mt . Using (2.1) in H lnt we 
obtain terms with both photon creation (a*) and atomic raising (Oi<j) operators (or photon destruction and 
atomic lowering {a, Oj>i}), or two photon creation (destruction) operator terms in Hq . Such high frequency 
terms lead to vacuum fluctuations and higher order processes. The rotating wave approximation sets these 
processes to zero. We also set to zero terms of the form tpi < jOk<i (or ^>] < _ 7 -Oj>fc) for ^ (k,l) since 

they give no contribution to lowest order in perturbation theory. So we only retain those terms which pair 
creation/lowering or destruction/raising operators (and creation/destruction operators in Hq) and connect 
states with w equal energy. 

Combining these two approximations we get 

#int = ~\\ (Ix^^tujj^Otja - ^e-^Ojie^^), (2.7) 

J i<j 

Ht = -2i f dx YjUiAj^ij- (2-8) 
J i<j 

The free field hamiltonian follows from the field action, which using 

\k e \ 2 « => \d x tpij\ < u>ij\if>ij\, \dti>ij\ < Uij\i>ij\, (2-9) 

approximates to 

J dxdt^({d t c}>) 2 - {d x 4>f) w 2i J dxdt ^w^tft + BtWij- (2-10) 

(In (2.8) we have dropped the quadratic term that appears in (1.21). For electric fields small compared 
with e/a§, this term is negligible in relation to ifmt-) The phases e ±laij and e ±tuJijX can be absorbed into 
{ipij , iplj } an d \Pi3i®ii\i<j respectively without changing the commutation relations. The time dependent 
phase e~ Vjiiit ( e + lw «*) cancels the time dependence of O i<: j (Oj > i) coming from the free atomic hamiltonian. 
Thus we can set i?Q tom to zero and consider the model defined by the complete hamiltonian 

H = -2ifdxJ2 Wij^d^a - l - I dxJ2 UijPij (^Oij - $.Oji) . (2.11) 
Finally we can rescale Vij and tpjj as 

which gives for the commutator (2.6) 

[ij l3 (x,t),ijl(x',t)] =5 lk S 3 i6(x-x') 7 (2.13) 
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and denning 



<i = TTlxPiiV^ = \l w ™ = w mm+i (1 < m < r), (2.14) 



2V2 

the hamiltonian (2.11) becomes 

H = —i dx^2 Tpijdxipij ~i dx^2 %ij (^ijOij ~ i>ijOji) ■ ( 2 - 15 ) 

The first interaction term, tpijOij, causes an atomic transition from a lower energy state |j) to a higher 
energy state |i), along with the absorption of a photon of energy w u^. The second term, ip\-Oji, causes a 
transition from a higher energy state |i) to a lower energy state |j), along with the creation of a photon of 
energy w Wy. 

From H we can obtain the Heisenberg operator equations of motion. Explicitly we find 



d t O k< i = - ^ uJijipijOkj + ^2 uj lk ipikOu + ^2 ^aWiPkt 

j>l i<k *<I 

i*k 

l-l 

~ ZkjipljOji + LOkliit J2 Hm ' (2.16a) 



i> k m=k 

3*1 



d t O l>k = d t O^ 



k<l 



^2 Ulj^lftjk + ^2 Uik^lkOu + ^2 ^illpilOik 
j>l i<k i<[ 

l-l 



- ^ujkj^kjOij + uj kl tp kl ^2 Km, (2.16b) 

j> k m=k 

3*1 

j>m+l 

- ^2 Uim+l(lpim+lO im+1 + lf>} m+1 O m+li ) + ^2 Uim(lpimO im + ^\ m O m i) 

+ J2 U m +lj{ll>m+ljO m +lj+i>i l +lj jm+l) (2.16c) 
j>m+l 

(d t + d x )ip k< i = uJ k iOi k (2.16d) 

{dt + d x )^ k<l = w kl O kl . (2.16e) 

Each term in the equations of motion for the atomic operators has a simple physical interpretation. For 
example, consider the raising operator O k< \. The first summation in (2.16a) is a sum over all atomic and 
photon operator pairs, where the photon field destroys a photon of energy w;<j, and the atomic (raising) 
operator causes a transition from the lower state |j) to the higher state |k). The change in energy of the 
system corresponding to the atomic/field pair being uo k i (as is for every other term in (2.16a)). If we think of 
a field destroying (creating) a photon of energy Wj<j as "connecting" atomic states |i) to |j) (|j) to |i)), along 
with Oi<j (Ojyi) connecting |j) to |i) (|i) to |j)) and H m connecting |m) to |m) and |m+ 1) to |m+ 1), then 
(2.16a) is, aside from c-number factors, a sum over all atomic/field pairs connecting |1) to |k) through some 
intermediate state, i.e., |1) — ► |j) — > |k). A similar interpretation follows for Oi >k and Ti. m - 
These equations of motion have a zero-curvature representation 

[d t + A u d x +A x ] =0, (2.17a) 

=> 9 t A x - d x A t = [A x , A t ], (2.17b) 
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where A x and A t are matrices of quantum operators given by 

A x = H Y, A ^nHP n + „ £ (O^ + O^) + £ 2 « H^S + ^ E^) 

i<j 

-T, E ( 2 - 18a ) 



Km,n<r 



= - E ^mn^^-E^'(-^«+^ii)> ( 218b ) 



1 < m , n < r i<j 



w ?<j = w i<Jfc + io 2 k<v (2.19) 



provided that the lu^s satisfy 

for any intermediate value of k, i < k < j. This is trivially satisfied by uju +1 , in which case there is no 
such k. Here A^ n is the inverse Cartan matrix, and {E?<j, Hf^} are matrices in any representation p of sl n 
satisfying (0.2). Requiring (2.17) to be valid for all values of the arbitrary spectral parameter fi, is equivalent 
to the equations of motion (2.16) (of course provided (2.19) is satisfied). The equivalence can be shown by 
making use of the commutation relations (0.2). From (2.16c) we can derive a more compact form for the 
equations of motion corresponding to the TL m basis operators 

r u 

d t n m = -j2 A mu E E s «( «^« +° ik ^ii)- ( 2 - 2 °) 

u— 1 u<l<nk— 1 

The constraint (2.19) arises in forming a zero-curvature representation for the field ({tpij, i>\j]) equations 
of motion, and reduces the number of free parameters to r = n — 1, these being {uj m }i<m<r- The definition 
(2.14) shows that the constraint is equivalent to the requirement that the spontaneous decay rate (for a 
single atom) from |i) to |j), Tfp be equal to the sum of the decay rates r* fc and T s k j for any intermediate 
state |k) 

T"<j — Ti<fc + ^i<j- (2.21) 

The equations of motion for the atomic operators have the zero-curvature representation (2.17) independent 
of the constraint. 

A zero-curvature representation implies that the system is integrable. Thus the model (2.15) can now 
be solved by the Quantum Inverse Scattering Method. 
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